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Abstract 

The production of the excited charm mesons Di(2420) and 1^2(2460) in ep colli- 
sions has been measured with the ZEUS detector at HERA using an integrated 
luminosity of 373 pb"^. The masses of the neutral and charged states, the widths 
of the neutral states, and the helicity parameter of Z}i(2420)'^ were determined 
and compared with other measurements and with theoretical expectations. The 
measured helicity parameter of the allows for some mixing of 5*- and D- 
waves in its decay to D*^7r^. The result is also consistent with a pure D-wave 
decay. Ratios of branching fractions of the two decay modes of the Z}2(2460)° 
and D*(2460)± states were measured and compared with previous measurements. 
The fractions of charm quarks hadronising into Di and D2 were measured and 
are consistent with those obtained in e~^e~ annihilations. 
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1 Introduction 



The production of the well-estabhshed ground-state charm mesons D and D* has been 
extensively studied in ep coUisions at HERA. The large charm production cross section at 
HERA makes it possible to also investigate the excited charm-meson states. In a previous 
ZEUS analysis [1], with an integrated luminosity of 126 pb~^, the orbitally excited states 
Di(2420)° with = 1+ and D;(2460)° with = 2+ were studied in the decay modes^ 
L'i(2420)° D*(2010)+7r- and L)*(2460)° D*(2010)+7r-, D+n-. The width of the 
Di was found to be significantly above the 2008 world-average value [2]. A study of the 
helicity angular distribution of the Z)i(2420)° gave results that were consistent with some 
5*- wave admixture in the decay — )■ D*^7r~, contrary to Heavy Quark Effective Theory 
(HQET) predictions [3,4] and to previous experimental results [5] which had yielded a 
pure D-wave decay in this channel. 

In this paper the analysis was repeated with an independent data sample of higher in- 
tegrated luminosity. In addition the production of the charged excited charm mesons 
Z)i(2420)+ and D2(2460)+ was studied for the first time at HERA in the decay modes 
^1(2420)+ ^ D*(2007)°7r+ and D*(2460)+ -> L'*(2007)°7r+, D°7r+. For both the neutral 
and charged excited charm mesons the study also includes a measurement of fragmenta- 
tion fractions and ratios of the D2 branching fractions. 

The analysis was performed using data taken from 2003 to 2007, when HERA collided 
electrons or positrons at 27.5 GeV with protons at 920 GeV. The data correspond to an 
integrated luminosity of 373 pb~^. The upgraded ZEUS detector included a microvertex 
detector, allowing the measurement of the decay vertex of charm mesons. In particular, 
the signal-to-background ratio was significantly improved for the meson, which has 
the highest lifetime among the charm hadrons. 

To maximise the statistics, both photoproduction and deep inelastic scattering events 
were used in this analysis. Events produced in the photoproduction regime contributed 
70 — 80% of the selected charm-meson samples. 

2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [6]. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central 
tracking detector (CTD) [7] and the microvertex detector (MVD) [8]. These components 

^ The corresponding anti-partiele decays were also measured. Hereafter, charge conjugation is imphed. 
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operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The 
CTD consisted of 72 cylindrical drift- chamber layers, organised in nine superlayers cover- 
ing the polar-angle^ region 15° < 9 < 164°. The MVD silicon tracker consisted of a barrel 
(BMVD) and a forward (FMVD) section. The BMVD contained three layers and pro- 
vided polar-angle coverage for tracks from 30° to 150°. The four- layer FMVD extended 
the polar- angle coverage in the forward region to 7°. After alignment, the single-hit 
resolution of the MVD was 24 |j.m. The transverse distance of closest approach (DCA) 
of tracks to the nominal vertex in the X-Y plane was measured to have a resolution, 
averaged over the azimuthal angle, of (46 © 122/pr) |^m, with px in GeV. For CTD- 
MVD tracks that pass through all nine CTD superlayers, the momentum resolution was 
(^{Pt)/pt = 0.0029pT © 0.0081 © 0.0012/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (GAL) [9] consisted of three parts: 
the forward (FGAL), the barrel (BGAL) and the rear (RGAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMG) and either one (in RGAL) or two (in BGAL and FGAL) hadronic sections (HAG). 
The smallest subdivision of the calorimeter was called a cell. The GAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/v^ for electrons and 
a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — )■ e^yp by a detec- 
tor which consisted of an independent lead-scintillator calorimeter [10] and a magnetic 
spectrometer [11] system. 

3 Event simulation 

Monte Garlo (MG) samples of charm and beauty events were produced with the Pythia 
6.221 [12] and the Rapgap 3.000 [13] event generators. The generation included direct 
photon processes, in which the photon couples directly to a parton in the proton, and 
resolved photon processes, where the photon acts as a source of partons, one of which 
participates in the hard scattering process. The GTEQ5L [14] and the GRV LO [15] 
parametrisations were used for the proton and photon parton density functions, respec- 
tively. The charm- and beauty-quark masses were set to 1.5 GeV and 4.75 GeV, respec- 
tively. The masses and widths for charm mesons were set to the latest PDG [16] values. 



^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
nominal proton beam direction, referred to as the "forward direction" , and the X axis pointing left 
towards the centre of HERA. The coordinate origin is at the centre of the CTD. The pseudorapidity 
is defined as rj — — In (tan |), where the polar angle, 9, is measured with respect to the Z axis. 
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Events for all processes were generated in proportion to the respective MC cross sections. 
The Lund string model was used for hadronisation in Pythia and Rapgap. The Bowler 
modification [17] of the Lund symmetric fragmentation function [18] was used for the 
charm- and beauty-quark fragmentation. 

The Pythia and Rapgap generators were tuned to describe the photoproduction and 
the deep inelastic scattering regimes, respectively [1]. Subsequently, the Pythia events, 
generated with < 1.5 GeV^, were combined with the Rapgap events, generated with 
> 1.5 GeV^, where is the exchanged-photon virtuality. 

The generated events were passed through a full simulation of the detector using Geant 
3.13 [19] and processed with the same reconstruction program as used for the data. 

4 Event selection and reconstruction of ground-state 
charm mesons 

The ZEUS trigger chain had three levels [6,20,21]. The first- and second-level trigger used 
GAL and GTD data to select ep collisions and to reject beam-gas events. At the third- 
level trigger, the full event information was available. All relevant trigger chains were 
used for the data. Triggers that required the presence of a reconstructed D*^ — t- D^tt^ — )■ 
[K~7r^)7i^ or (A'^7r+7r^7r+)7r+, — )■ /i^7r+7r+ or — )■ K^n^ candidate constituted 
a major fraction of the selected events. However, events missed by these triggers but 
selected with other trigger branches were also used in the analysis. Applying, in the MG, 
either no trigger selection cuts or requiring at least one trigger chain to be passed did not 
affect the final measurements. 

To ensure high purity in the event sample, the Z position of the primary vertex, recon- 
structed from GTD and MVD tracks, had to be within |Zvtx| < 30 cm. All charm mesons 
were reconstructed with tracks measured in the GTD and MVD. All tracks were required 
to have a transverse momentum, px, above 0.1 GeV, to start not further out than the 
first GTD superlayer and to reach at least the third superlayer. The tracks were assigned 
either to the reconstructed primary vertex or to a secondary decay vertex associated with 
the weak decay of a charm meson, or D^. To ensure the use of well reconstructed 
MVD tracks, all tracks associated with the secondary vertex were required to have at 
least two BMVD measurements in the X-Y plane and two in the Z direction. 

The decay-length significance is a powerful tool for rejection of combinatorial background. 
It is defined as S = l/ai, where the decay length / is the distance in the transverse 
plane between the production point and the decay vertex of a candidate charm meson 
projected on its momentum direction and ai is the uncertainty of this quantity [22]. The 
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quantity S is positive when the angle between the particle momenta and the direction from 
primary to secondary vertex is less than 7r/2; it is negative otherwise. The S distribution 
is asymmetric around zero, with a stronger positive contribution coming mostly from 
the charm mesons. The contributions to negative S values are due to background and 
resolution effects. 

The combinatorial background was suppressed by selecting events above a minimum value 
of the ratio pt{D) / E^^^° , where D denotes D*~^, or and E^^^^° is the transverse 
energy measured using all CAL cells outside a cone of 10° around the forward direction. 
In addition, to reduce background, the dE/dx values measured in the CTD of track 
candidates originating from the D mesons were used. The parametrisation of the dE/dx 
expectation values and the probabilities Ik and /^r of the kaon and pion hypotheses, 
respectively, were obtained as described in previous analyses [23,24]. The cuts Ik > 0.03 
and It, > 0.01 were applied. 

4.1 Z)*+ reconstruction 

D*~^ mesons were identified via the decay modes D*~^ D^'^t ~^ {K~7t~^)7t'^ and D*~^ — )> 
D'^TT^ — (Ji ~7r+7r~7r+)7r^, where tTs is a low-momentum ("soft") pion due to the small 
mass difference between D*~^ and D^. Tracks were combined to form candidates by 
calculating the invariant-mass combinations M^Kn) or M[K'K'kt[) with total charge zero. 
D*^ candidates were formed by adding a soft pion, tt^, with opposite charge to that of the 
kaon. Combinatorial background was reduced by applying cuts as detailed in Table 1. 

The mass differences AM = M{KTnT,) - M{K-k) and AM = M{K7nnnis) - M{K7nni) 
were calculated for the D*~^ candidates that passed the cuts of Table 1. Figure 1 shows 
the AM distributions for these D*~^ candidates. Clean peaks are seen at the nominal 
value of M{D*+) - M{D^) [16]. 

The AM distributions were fitted to a sum of a background function and a modified 
Gaussian function [1] . The fit yielded D*^ signals of 64988 ± 430 candidates for — > 
i^TT and 24441 ±310 candidates for — )■ Kti-kti. The fitted mass differences were 
145.400 ± 0.003 MeV and 145.420 ± 0.003 MeV respectively, in agreement with the PDG 
average value [16]. Only D*'^ candidates with 0.144 < AM < 0.147 GeV were used for 
the excited charm mesons analysis. 

4.2 reconstruction 

mesons were reconstructed from the decay — t- A'^tt+tt^ with looser kinematic 
cuts than in the previous analysis [1], made possible by the cleaner identification with 
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the MVD. For each event, track pairs with equal charge and pion mass assignment were 
combined with a track with opposite charge with a kaon mass assignment to form a 
candidate. These tracks were refitted to a common decay vertex, and the invariant 
mass, M{K7r7r), was calculated. The K and vr tracks were required to have transverse 
momentum pff > 0.5 GeV and p'^ > 0.35 GeV and the distance of closest approach 
between each pair of the three tracks was required to be less than 0.3 cm. To suppress 
combinatorial background, the following cuts were applied: 

• cos 6'*(i^) > —0.75, where 6*{K) is the angle between the kaon in the Kim rest frame 
and the Kinr line of flight in the laboratory frame; 

• the of the fit of the decay vertex was less than 10; 

• the decay-length significance, S{D~^), was greater than 3. 

Background from D*~^ decays was removed by requiring M^Kttti) — M{Kti) > 0.15 GeV. 
Background from — )■ 07r, — )■ K^K~ was suppressed by requiring that the invari- 
ant mass of any two decay candidate tracks with opposite charge should be outside 
±8 MeV around the nominal mass when the kaon mass was assigned to both tracks. 

candidates in the kinematic range pt{D^) > 2.8 GeV and \ri{D~^)\ < 1.6 were kept 
for further analysis. 

Figure 2 (a) shows the M{K~7r^7r^) distribution for candidates after the cuts. A clear 
signal is seen at the nominal value of the mass [16]. The mass distribution was fitted 
to a sum of a modified Gaussian function and a polynomial background. The fit yielded 
a signal of 39283 ±452 events and a mass of 1869.1 ± 0.1 MeV, in agreement with 
the PDG average value [16]. Only candidates with 1.85 < M^Ktitt) < 1.89 GeV were 
used for the excited charm mesons analysis. 

4.3 reconstruction 

mesons were reconstructed from the decay — j- A'"7r+. For each event, two tracks 
with opposite charge and K and vr mass assignments, respectively, were combined to form 
a candidate. These tracks were refitted to a common decay vertex, and the invariant 
mass, M{K7i), was calculated. Both tracks were required to have transverse momentum 
p!f > 0.5 GeV and pj. > 0.7 GeV and the distance of closest approach between these 
tracks was required to be less than 0.1 cm. To suppress combinatorial background, the 
following cuts were applied: 

• I cos 6'*(ii') I < 0.85, where 6*{K) is the angle between the kaon in the Kn rest frame 
and the Kn line of fiight in the laboratory frame; 

• the of the decay vertex was less than 20; 
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• the decay-length significance, 5'(D°), was bigger than 0. 

D'^ candidates in the kinematic range pt{D^) > 2.6 GeV and \r]{D^)\ < 1.6 were kept for 
further analysis. 

Figure 2 (b) shows the M{K~7r^) distribution for candidates after the cuts. A clear 
signal is seen at the nominal value of the mass [16]. The mass distribution was 
fitted to a sum of a modified Gaussian function, a broad modified Gaussian representing 
the reflection produced by mesons with the wrong (opposite) kaon and pion mass 
assignment and a polynomial background. For the reflection, the shape parameters of the 
broad modified Gaussian were obtained from a study of the MC signal sample and the 
normalisation (integral) was set equal to that of the other modified Gaussian. The fit 
yielded a signal of 145740 ± 2944 events and a D° mass of 1864.1 ± 0.1 MeV which 
is 0.8 MeV lower than the PDG average value [16]. This deviation does not affect any of 
the results of the excited charm mesons. Only candidates with 1.845 < M^Kttit) < 
1.885 GeV were used for the excited charm mesons analysis. 



5 Di and reconstruction 

5.1 Reconstruction of the Di(2420)^ and D*{2460f mesons 

The Di{2A20)^ and 1^2(2460)° mesons were reconstructed in the decay mode D*^7t^ by 
combining each D*'^ candidate with an additional track, assumed to be a pion (tTq), with 
a charge opposite to that of the D*. Combinatorial background was reduced by applying 
the following cuts: 

• PriTTa) > 0.15 GeV; 

• ri{7ia) < 1.1 ; 

• pT{D*+7ia)/E'^>^'^° > 0.25 (0.30) for the Kn (D° Kmni) channel; 

• cos9*{D*+) < 0.9, where 9*{D*'^) is the angle between the Z)*+ in the D*^-Ka rest 
frame and the D*'^Tia line of flight in the laboratory frame; 

• the cut > 0.01 was applied for tTq. 

For each excited charm-meson candidate, the "extended" mass difference, AM*^^* = 
MiKTriTsiTa) - M{KTnTs) or AM^^* = M{KmnTnsiTa) - M{K'inxTxHs), was calculated. 
Figure 3 (a) shows the invariant mass M(D*+7ra) = AM'=''* + M(D;+q), where M(L';+q) 
is the nominal D*^ mass [16]. A clear signal in the D\/ D"^ mass region is seen. 

The 1^2(2460)° was also reconstructed in the decay mode ^2(2460)" — t- D^ti^ by com- 
bining each candidate with an additional track, assumed to be a pion tTq, with a 
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charge opposite to that of the D^. Combinatorial background was reduced by applying 
the following cuts: 

• prin,) > 0.3 GeV; 

• r]{TTa) < 1.5; 

. pT{D^na)/E'>'^° > 0.35; 

• cos6*{D^) < 0.8, where 9*{D^) is the angle between the in the D'^Ha rest frame 
and the D^iTa line of flight in the laboratory frame; 

• the cut > 0.01 was applied for tTq. 

The D2(2460)° — D'^rc^ decay mode was reconstructed by calculating the "extended" 
mass difference AM^^* = M{K7i7r7Ta) — M{K7nr). Figure 3 (b) shows the invariant mass 
M{D+7ra) = AM^^* + M(L)^g), where M(L)^g) is the nominal D+ mass [16]. A 
clear signal is seen. No indication of the — > D^ir^ decay is seen, as expected 
from angular momentum and parity conservation for a = 1"*" state. The various 
contributions to the mass spectrum will be discussed below. 

5.2 Reconstruction of the £>i(2420)+ and D*{2460)+ mesons 

The charged excited meson Di(2420)"'" has been seen [16] in the decay modes D*^7t~^ 
and D^n^Ti^ and the charged excited meson 1^2(2460)"^ has been seen [16] in the decay 
modes D*^it^ and D^7t~^. A search for Df and 1^2^ signals was performed in the mass 
distribution M{D^7r'^). For the a possible D^vr"*" signal can arise only via a feed- 
down contribution (see Section 6). Each candidate was combined with an additional 
track, assumed to be a pion (vTa), with either positive or negative charge. Combinatorial 
background was reduced by applying the following cuts: 

• PtM > 0.35 GeV; 

• vi^^a) < 1-6; 

. pTiD%a)/El>''" > 0.3; 

• cos^*(-D'^) < 0.85, where 9*{D^) is the angle between the in the -D'^tTq rest frame 
and the D^iTa line of flight in the laboratory frame; 

• the cut > 0.01 was applied for tTq. 

For each excited charm-meson candidate, the "extended" mass difference AM°^* = 
M[KmTa) — M{K'k) was calculated. Figure 4 shows the invariant mass M{D^Tra) = 
AM^^* + M(L>^j3g)' where M{D^j^q) is the nominal mass [16]. A clear signal of 
Z^2^ —7- D^Ti^ is seen. An enhancement above background is also seen at the mass region 
around 2.3 GeV. The various contributions to the mass spectrum will be discussed below. 
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6 Mass, width and helicity parameters of Di and X)J 



A significant enhancement above background is seen in the D^7r~^ mass distribution (Fig. 4) 
around 2.3 GeV. A small excess of events is also seen in the same mass region in the D^tt~ 
mass distribution (Fig. 3(b)). 

The origin of these structures in both spectra is similar. They originate from the de- 
cay chains D^^,Df D*+7i-, D*+ D+7i° and Dt,D;+ D*°7i+,D*^ D°7i° or 
£)*o _s, £)0^_ rjj-^g 71*^/7 are not seen in the tracking detectors; thus, the reconstruction 
is incomplete. However, since the available phase space in the D* — )■ Dir^ decay is small 
and D is much heavier than tt", the energy and momentum of D are close to those of D*. 
Consequently, the enhancements in the MlD^iTa) (Fig. 3(b)) and M{D^TTa) (Fig. 4) distri- 
butions are feed-downs of the excited charm mesons Di,D2, shifted down approximately 
by the value of the 7r° mass, as verified by MC simulations. 

6.1 Fitting procedure for and D*^ 

To distinguish between D\^D^ D*^7r~, their helicity angular distributions were used. 
These can be parametrised as dN/ d cos a (x.l + h cos^ a, where a is the angle between the 
tTq and vr^ momenta in the D*^ rest frame and h is the helicity parameter, predicted [3,4] 
to be /i = 3 for D\ and h = —1 for D"^. Figure 5 shows the M{D*^TTa) distribution in 
four helicity bins. As expected from the above h values, the D\ contribution increases 
with I cosa| and dominates for | cosa| > 0.75, where the D"^ contribution is negligible. 

A fit was performed using simultaneously the M{D~^'jTa) distribution (Fig. 3(b)) and 
the M{D*'^7Ta) distributions in four helicity bins (Fig. 5). The background was de- 
scribed by four parameters a,b,c,d, separately for M{D*~^'7ia) and M{D~^Ha), as B{x) = 
ax^exp(— cx — dx"^), where x = AM^^^ — M^+. Each resonance was fitted to a relativistic 
Z)-wave Breit-Wigner (BW) function [1] convoluted with a Gaussian resolution function 
with a width fixed to the corresponding MC prediction. Yields of the three signals, the 
and Df^ masses and widths and the helicity parameter, h{D^), were free parameters 
of the fit while h{D'^) was fixed to the theoretical prediction [3,4], h{D'^) = — 1. Another 
free fit parameter was the contribution of the Df^ feed-downs to the M{D'^'Ka) dis- 
tribution (see Appendix). The total normalisation of the sum of the feed-down processes 
from and decays was fitted relative to the direct signal peak yield from decay. 
The relative yields of the two feed-down contributions were taken to be equal to those for 
the direct signals in the D*~^7t~ decay channel. 

The wide excited charm states [16] -Di(2430)° and -00(2400)" are expected to contribute 
to the M{D*~^7ra) and M{D~^7ra) distributions, respectively. Even though these states are 
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hardly distinguishable from background due to their large width, they were included in 
the simultaneous fit with shapes described as relativistic S'-wave BW functions [1] . Their 
masses and widths were set to the PDG values [16]. The yield of the 1^1(2430)° was set 
to that of the narrow Di(2420)° meson since both have the same spin-parity = 1+. 
The ratio of Dq{2400)^ to the narrow state D2(2460)° was a free parameter in the fit. 

The results of the simultaneous fit are given in Table 2 and shown in Figs. 3 and 5. 
Systematic uncertainties are discussed in Section 8. All results from the new analysis 
(HERA II) are consistent with those from the previous ZEUS publication [1] (HERA I). 
The masses of both and -Dg" are consistent with the PDG values [16] and with a recent 
BABAR measurement [25]. The width, T{D^) = 38.8 ± 5.0(stat.)+5:4(syst.) MeV, is 
also consistent with the PDG value [16] of 27.1 ±2.7 MeV, and is in good agreement with 
the BABAR measurement [25] of 31.4 ± 0.5 ± 1.3 MeV. The Df width, T{Df) = 46.6 ± 
8.1(stat.)+|^(syst.) MeV, is consistent with the PDG value [16] of 49.0 ± 1.4 MeV, and 
with the BABAR measurement of50.5±0.6±0.7 MeV. 

The helicity parameter, h{D'^) = 7.8+2;7(stat.)+i:^(syst.), is consistent with the BABAR 
value of h{Dl) = 5.72 ± 0.25 and somewhat above the theoretical prediction of /i = 3 and 
measurements by CLEO [26] with h{D^) = 2.74^1^1- The simultaneous fit with h{D^) 
fixed to the theoretical prediction, h{D\) = 3, yielded masses and widths of -D^^ and Z)° 
that are somewhat away from the PDG values [16]. Repeating the simultaneous fit with 
h{D^) as a free parameter yielded similar results for all other free parameters with some- 
what larger errors and with h{D^) =—1.16 ±0.35, in good agreement with the theoretical 
prediction oi h = —1. 

The helicity angular distribution for a J'^ = 1"^ state with a mixture of D- and S'-wave is 
dN 



(icosa 



ocr±(l-r)(l±3 cos^ a) /2 ± v^2r(l - r) cos 0(1-3 cos^ a) , (1) 



where r = r5'/(r5' ± F/j), F5(F£)) is the S{D)-wave partial width and is relative phase 
between the two amplitudes. The relation between h, r and is given by 

(3-ft)/(3 + /.)-r 

The range of the measured h{D\) restricted to one standard deviation is shown in Fig. 6 
in a plot of cos0 versus r. This range is consistent with the BABAR measurement [25]. 
The range restricted by CLEO [26] is outside the range of this measurement and that 
of BABAR. A similar measurement by the BELLE collaboration [5] is consistent with a 
pure D-wave, i.e. F5/(F5 ± Vjj) = 0. 

In a recent paper [25] the BABAR Collaboration searched for excited D meson states 
in e~^e~ cc ^ D'^*\ + X with very large statistics. In addition to the D\ and 
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resonances, they saw two new structures near 2.6 GeV in the D^7r~ and D*~^7r~ mass 
distributions, D(2550)° and -D*(2600)°, and interpreted them as being radial excitations 
of the well-known and Z)*°, respectively. A small enhancement of events above the 
solid curve in the region near 2.6 GeV is seen in the M{D*~^7t~) distribution (Figs. 3(a), 5). 
Adding the new BABAR states to the fit gave insignificant yields of the states and did 
not significantly change the results of the other fit parameters. 

6.2 Fitting procedure for and 

To extract the Df and D^^ masses and yields, a minimal fit w^-s performed using 
the M(D°7r„) distribution (Fig. 4). Both resonances were fitted to relativistic D-wave 
Breit-Wigner (BW) functions [1] convoluted with a Gaussian resolution function with a 
width fixed to the corresponding MC prediction. Yields of the -D2^ — )■ D^7r~^ and the 
two feed-downs DfjD'^'^ — )■ D*^7r~^ (see Appendix) and the Df and masses were 
free parameters of the fit. The Df and -Dg^ widths were fixed to the PDG values [16] 
and the Df and Dp~ helicities were fixed to the theoretical prediction [3,4], h{Df) = 3 
and h{D2'^) = —1. The background was parametrised with four parameters a,b,c,d as 
B{x) = ax^ exp{—cx — dx^), where x = AM°^* — M^+. 

The results of the fit (yields and masses) are given in Table 3 and shown in Fig. 4. The 
masses of Df and are consistent with the PDG values [16]. The mass is also 
consistent with the BABAR measurement [25]. 

7 D2 decay branching ratios and D1/D2 fragmenta- 
tion fractions 

7.1 The neutral excited mesons 

The branching ratio for and the fragmentation fractions for and were mea- 
sured using the channels D+tt' and D?, D^^ D*+7r- with D*+ D°it+ 
(i^"7r+)7r+. The numbers of reconstructed D^,Df D*^7r~ and D2° ~^ D^n^ decays 
were divided by the numbers of reconstructed D*"^ and mesons, yielding the fractions 
of D*~^ and mesons originating from the and 1^2° decays. To correct the measured 
fractions for detector effects, ratios of acceptances were calculated using the MC simu- 
lation for the D^, Df D*+7r" and Df D^n^ states to the inclusive and 
acceptances, respectively. 

Beauty production at HERA is smaller than charm production by two orders of magni- 
tude. A subtraction of the 6-quark relative contribution in a previous ZEUS analysis [1] 
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changed the relative acceptances by less than 1.5% of their values. Consequently, no such 
subtraction was performed in this analysis and the MC simulation included the beauty 
production processes. A variation of this contribution was considered for the systematics 
(Section 8). 

The fractions, J-", of D*'^ mesons originating from D\ and D'^ decays were calculated in 
the kinematic range \rj{D*^)\ < 1.6 and pt{D*~^) > 1.5 GeV for the D*~^ decay and the 
fraction of mesons originating from decays was calculated in the kinematic range 
Pt{D+) > 2.8 GeV and \viD+)\ < 1.6. 

The fractions measured in the restricted pt{D*^ , D^) and ri{D*^ , D^) kinematic ranges 
were extrapolated to the fractions in the full kinematic phase space using the Bowler 
modification [17] of the Lund symmetric fragmentation function [18] as implemented in 
Pythia [27]. Applying the estimated extrapolation factors, ~ 1.12 for d1^d*+tt- /d*+-i 
~ 1.16 for J^Df^D'+^-/D'+ and ~ 1.34 for d*^o^d+^- /d+, gives 

nT-.D*+.-/D*+ = 8-5 ± 1.4(stat.)l}:^(syst.) %, (3) 
n1o_z,.+.-/z,.+ = 4.7 ± 1.3(stat.)lJJ(syst.) %, (4) 

•^D;o%m.-/D+ = 6.7 ± 2.4(stat.)+}i(syst.) %. (5) 

In the full kinematic phase space, the extrapolated fractions of D*^ originating from 
D\ and D'^ and of originating from D"^ can be expressed [1] in terms of the rates 
of c-quarks hadronising to a given charm meson ("fragmentation fractions"), /(c — )■ -D?), 
f^c-^D^^), f{c^D*~^) and f{c^D^) and the corresponding branching fractions 

^Dl->D*+Tz-^ ^Df~^D*+n- and Bjj*0^jj+^-. 

From the expressions used in a previous ZEUS publication [1], the fragmentation fractions 
f{c^ D^) and f{c^ D'^) and the ratio of the two branching fractions for the meson 
can be shown to be: 

/(c ^ /^?) = ^f"""--^"" /(c ^ /^-), (6) 



(8) 
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The /(c — !■ D*~^) and /(c — D"*") values used were obtained as a combination of data from 
HERA and e+e^ colliders [28]: 

/(c ^ D*+) = 22.87 ± 0.56(stat. © systO^aselbr.) %, 

/(c ^ D+) = 22.56 ± 0.77(stat. © syst.) ± 1.00(br.) %. 

where the third uncertainties are due to the branching-ratio uncertainties. 

Taking into account the correlations in the simultaneous fit performed to obtain the values 
in Eqs. (4) and (5) yields 

Df^D+n- ^ i.4±o.3(stat.)±0.3(syst.), 

in good agreement with the PDG world-average value 1.56 ± 0.16 [16]. Theoretical mod- 
els [29-31] predict the ratio to be in the range from 1.5 to 3. 

Neglecting the contributions of the non-dominant decay mode — )■ D^ti^ti^ [16] and 
assuming isospin conservation, for which 

^D'lWD' + TT- = 2/3, Bjj,0^j;), + ^- + Bj;)*0_,jj+^- = 2/3, 

and using Eqs. (6) and (7) yields 

/(c ^ L>?) = 2.9 ± 0.5(stat.) ± 0.5(syst.) %, 

/(c ^ Df) = 3.9 ± 0.9(stat.)+°:^(syst.) %. 

The measured fragmentation fractions were found to be consistent with those obtained in 
e~^e~ annihilations [32]. The sum of the two fragmentation fractions, 

/(c ^ + /(c ^ D;^) = 6.8 ± 1.0(stat.)+[J:^(syst.) %, 

agrees with the prediction of the tunneling model of 8.5% [33]. 
Assuming uncorrelated errors, the ratio 

/(c ^ Dl)/f{c Df) = 0.8 ± 0.2(stat.) ± 0.2(syst.) 

is in good agreement with the simple spin-counting prediction of 3/5. 
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7.2 The charged excited mesons 



The branching ratio for 1^2^ and the fragmentation fractions for and were mea- 
sured using the channels D^ti+ and , D;+ D*°7r+ with D*° D^ny-f, 
where the 7r°/7 are not measured directly. Since D*^ decays always to [16], the 
number of D*^ and originating from / D^^^ are identical. The number of recon- 
structed Dt/D*+ D*^n+; D^n'^/-f and D°7r+ decays were thus divided 
by the total number of reconstructed mesons, yielding the fractions of mesons 
originating from Df /D2~^ decays. Detector effects were corrected as described in Sec- 
tion 7.1. The above fractions were calculated in the kinematic range pt{D^) > 2.6 GeV 
and \ri{D^)\ < 1.6 and extrapolated to the fractions in the full kinematic phase space 
as for the and (Section 7.1). Applying the extrapolation factors, ~ 1.28 for 
Df D*°7r+, ~ 1.18 for D*^7i+ and ~ 1.35 for D^n+ gives 

^Ti^D*^.VD^ = 5-4 ± 2.1(stat.)12;3(syst.) %, (9) 

•^d1+^d*o.+/z,o = 1-8 ± 0.9(stat.)+°J(syst.) %, (10) 

-^D*+^DO.+/Do = 2.0 ± 0.5(stat.)l°:^(syst.) %. (11) 
The fractions of D*^/D^ mesons originating from / decays can be expressed as 

™tr ^ N{Dt ^ D*'7r+) _ f{c ^ Dt) ^ 

™tr ^ iV(j;+ ^ D-°7r+) _ /(c ^ D;+) ^ 

where N denotes the acceptance-corrected number of events. 

The ratio of the fragmentation fractions /(c — j- D*^) and /(c — )■ D°) can be expressed as 

/(c ^ D*°) iV(D*0) 



/(c^DO) iV(Z}0) 
Consequently, Eqs. (12) and (13) can be written as 
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N{D;+ D*^Ti+) _ f{c Dp 



yielding 



Neglecting the non-dominant decay mode D'^ — )■ D^tt^h [16], assuming isospin conser- 
vation, for which 

^D+->D*07r+ = 2/3, 'S^*+_^£)*0^+ + ^D*+^D'-^TT+ = 2/3, 

and using Eqs. (9-11) and the fragmentation fraction [28] 

/(c ^ D^) = 56.43 ± 1.51(stat. © syst.)^J:64(br.) %, 



gives 



/(c ^ D+) = 4.6 ± 1.8(stat.)l^;^(syst.) %, 
/(c -> = 3.2 ± 0.8(stat.)l|];^(syst.) %, 
/(c ^ D+) + /(c ^ = 7.8 ± 2.0(stat.)+2;0(syst.) %, 



/(c -> /^i+)//(c ^2*+) = 1.4 ± 0.7(stat.)^°i(syst.), 

in agreement with the fragmentation fractions of the neutral excited charm mesons (Sec- 
tion 7.1). 
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The ratio of the branching fractions of the two dominant decay modes of the D2 



D*+^DOn+ _ ^ o.4(stat.)lo:2(syst.), (15) 



^D*+->D'-On+ 



significantly improves on the accuracy of the PDG [16] value of 1.9 ± 1.1 ± 0.3. BABAR 
measured the ratio [34] ^ = 0.62 ± 0.03 ± 0.02, which depends on 



some assumptions and is not included in the PDG averages [16]. Using the value given 

in Eq.(15) yields a ratio 75 1^2^^^°-+ ^ 0.52+°-°^(stat.) ± 0.05(syst.), in good 

agreement with the BABAR result. 



8 Systematic uncertainties 



The systematic uncertainties were determined by appropriate variations of the analysis 
procedure, generally by the uncertainties in our knowledge of the variables considered, 
and repeating the calculation of the results. The following sources of uncertainty were 
considered: 

• {61} The stability of the fit results was checked by a variation of the selection cuts 
which are most sensitive to the ratio of signal and background in the data: 

- the cut on the minimal transverse momentum of the D*^, and candidates 
was varied by ±100 MeV; 

- the cut on the minimal transverse momentum of the extra pion in the excited D 
meson analysis was varied by ±10 MeV; 

- the selection cut on the cosine of angle between extra pions and charged (neutral) 
excited D meson candidates was changed by ±0.1 (±0.05); 

- the widths of the mass windows used for the selection of D*^ and candidates 
in the excited charm meson analyses were varied by ±5% for each pt dependent 
window (see Table 1), while for the candidates it was varied by ±12.5%. 

• {62} The CAL energy scale is known with ±2% uncertainty and was varied accordingly 
in the simulation. 

• {63} The uncertainties related to the fit procedure were obtained as follows: 

- the ranges for the signal fits were reduced on either side by 16 MeV for the D*~^7i 
and D'^Ti mass spectra and 24 MeV for the D^ir mass spectrum; 

- the background shape was changed to that used by BABAR (Eq. 1 in ref. [25]); 
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- the widths of the Gaussians used to parametrise the mass resolutions were changed 
by ±20%; 

- all the masses and widths of wide states were set free in the fit. Since with the 
present data alone these parameters are not determined well, the world-average 
values from PDG [16] were used as additional constraints. This was implemented 
by adding for each parameter P (width or mass) a term *-^^"p^™°] to the x^- 
function. Here Ppdg and o"(Ppdg) denote the parameter value and its uncertainty 
from PDG [16]; 

- the background functions in the four helicity intervals were allowed to have sepa- 
rate normalisations; 

- the helicity parameter of the Dg*^ meson in the fit was set free (Section 5.1). 

• {64} The uncertainties of M{D*'^)pdq, M(D°)pdg5 ^(-D'^)pdg were taken into ac- 
count. 

• {S^} The widths of Df and D^'^ were varied within their uncertainties taken from 
PDG [16]. 

• {Sq} The uncertainty of the beauty contamination was determined by varying the 
beauty fraction in the MC sample between and 200% of the reference amount. 

• {67} The extrapolation uncertainties were determined by varying relevant parameters 
of the Pythia simulation using the Bowler modification [17] of the Lund symmetric 
fragmentation function [18]. The following variations were performed: 

— the mass of the c quark was varied from its nominal value of 1.5 GeV by ±0.2 GeV; 

— the strangeness suppression factor was varied from its nominal value of 0.3 by 
±0.1; 

— the fraction of the lowest-mass charm mesons produced in a vector state was varied 
from its nominal value of 0.6 by ±0.1; 

— the Bowler fragmentation function parameter was varied from the predicted 
value 1 to 0.5; the a and b parameters of the Lund symmetric function were varied 
by ±20% around their default values [27]. 

A possible model dependence of the acceptance corrections was checked by reweighting the 
D-meson transverse momentum distribution in the MC to match the distribution observed 
in the data; no significant effect on any result was found. As a further cross check the 
selected pseudorapidity range of the extra pion, which is not the same for the different 
decay channels (see Section 5), was varied, and again no significant effect on any result 
was observed. The uncertainties of the fragmentation fractions /(c — ?■ D*^), f{c — )■ D^) 
and /(c — 7- D^) were included by adding in quadrature their statistical and systematic 
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uncertainties and the uncertainties originating from the branching-ratio uncertainties. 
The resulting uncertainty is included in 5-j. 

The contributions from all systematic uncertainties were calculated separately for posi- 
tive and negative variations and added in quadrature. The obtained values are listed in 
Tables 4-7. There is no single dominating source of systematic uncertainty. The total 
systematic uncertainties are comparable to the statistical errors. 

9 Summary 

The full HERA data taken from 2003 to 2007 with an integrated luminosity of 373 pb~^ 
has been used to study the production of excited charm mesons. Signals of Di(2420)° 
and 1^2(2460)'' were seen in the D*^ti~ decay mode and a clear Z}2(2460)° signal was seen 
in the D^t[~ decay mode. The measured and masses and widths are in good 
agreement with the latest PDG values. The measured helicity parameter allows for 
some S'-wave mixing in its decay to D*'^tt~ . The result is also consistent with a pure 
D-wave hypothesis. The helicity of Dj", when set free in the fit, is consistent with the 
HQET prediction, h = —1. 

A clear 1)2(2460)"'' signal is seen for the first time at HERA in the D^n'^ decay mode. 
Feed-downs of both resonances Z)i(2420)"'" and 1^2(2460)"'" in the decay mode D*^tt^ are 
seen in the expected mass region of M{D^7r'^) ^ 2.3 GeV. The measured Df and D2'^ 
masses are in good agreement with the PDG values and the D^^^^ mass is consistent with 
the BABAR measurement. 

The fractions of c-quarks hadronising into and ^ire consistent with those from 
the previous ZEUS publication and with e~^e~ annihilation results, in agreement with 
charm fragmentation universality. The fractions of c-quarks hadronising into Df and 
were measured for the first time and are consistent, respectively, with the fractions 
of the neutral charm excited states and Dg^- 

The ratios of the neutral and charged D2 branching ratios into Dn and -D*7r are consistent 
with the PDG values. 
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10 Appendix: Parametrisation of the feed-down con- 
tributions 

Let us consider the decay chain D12 — D*it, D* — t- Dir^ in the D* centre-of-mass 
system. Here D12 is a neutral (positively charged) excited charm meson Di or D2, D* is 
a positively charged (neutral) Z?*, tt is a negatively (positively) charged pion and _D is a 
positively charged (neutral) D (charge conjugation is implied). In this system Z)i.2 and tt 
in the initial decay and D and vr*^ in the subsequent decay are produced with back-to-back 
momenta. The momenta of particles in this system are: 



M2 - Ml, - MlV 



TT' 



where M is the Di 2 mass; 



1' Ml, - Ml + M%V 



The measured M{D7i) is given by 

Ml = M\D7r) = Ml + Ml + {P^ + Ml){P^ + M^) - 2PdP. cos a, 

where a is the helicity angle between 7r° and n. Using the equations above, can be 
parametrised as: 

Ml = M'^{1 - a) + b + gy^{M^ - di){M^ - rfs) cosa, (16) 



where 



a = (Mz).' + M^o' - Mn^)/i2MD*^), 
b = M,o2 - {Mn*' - M^^){Md^' + M^ - Mn')/{2Mn*') 

g = ^J {Md*^ + My - Mo^ - 4M^*'M,oV(2Mb.'), 
di = {Mn* + M^)^ 
d2 = {Md* - M^. 
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From Eq.(16), M is obtained as a function of Mm and a 

M = M{Mm,a). 

If the spectrum shape of M is 

^^/(-)' 

where TV is the number of candidates, then the Mm spectrum shape is 

r,..,.. NX dM 

f{M{Mm))- 



dMm ' ^'dMm 

Combining Eq.(16) with the normahsed hehcity angular distribution 

dN 1 + h cos^ a 
d{cosa) ^ 2{l + h/3) ' 

yields 

d'^N dM 1 + hcos"^ a 

f{M{Mm,a))- 



dMmd{cosa) ■> ^ ^ " dMm 2{l + h/?>) 
The fit uses the integral over cos a 

1 

dN XX dM l + /icos^a 



dMr, 



/^/,^/,, XX dM I + /1C0S a X /,«s 

/(M(M„..a))— ^^^^d(cosa). (17) 



Here f{M) is parametrised by a relativistic Breit-Wigner function as for the prompt 
signals. 

For the description of the D°7r spectrum, the D*^ D'^'j decay was also taken into 
account by replacing M^^o with My = in the equations above. For the description of the 
D~^7r spectrum, the contribution of the D*^ — )■ D'^'-f decay was neglected [16]. 
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Variable 


^0 ^ K-TI+ 




Pt{K) (GeV) 


> 0.45 


> 0.3 


PriTi) (GeV) 


> 0.45 


> 0.3 


priTTs) (GeV) 


> 0.1 


> 0.1 


PTiD*+) (GeV) 


> 1.5 


> 3 


\viD*^)\ 


< 1.6 


< 1.6 


Pt{D*+)/E'>^''° 


> 0.12 


> 0.18 


M{D^) (GeV) for 
Pt{D*+) < 3.25 GeV 


1.83 - 1.90 


1.84- 1.89 


M{D^) (GeV) for 
3.25 < Pt{D*+) < 5 GeV 


1.82 - 1.91 


1.84- 1.89 


M(L)O) (GeV) for 
5 < Pt{D*+) < 8 GeV 


1.81 - 1.92 


1.84- 1.89 


M{D^) (GeV) for 
Pt{D*+) > 8 GeV 


1.80 - 1.93 


1.84- 1.89 



Table 1: Cuts on D*^ — t- D^tx^ candidates for the decay channels — K vr 
and — )■ K'-k^ti^ti^ . 
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KFR A TT 


KFR A T 






2732 ± 285 


3110 ± 340 




i\[D2 — > U Tl) 


1798 ± 293 


870 ±170 




N{Df D+tt) 


521 ±88 {S{D+) > 3) 


690 ± 160 




M{D^), MeV 


2423.1 ± 1.5^?:^ 


2420.5 ±2.1 ±0.9 


2421.3 ±0.6 


r(DO), MeV 


38.8 ±5.0^^:^ 


53.2 ± 7.2+^1 


27.1 ± 2.7 




7 0+6.7+4.6 
' •°-2.7-1.8 


c- q+3.0+2.4 
'-'•^-1.7-1.0 




M{Df), MeV 


2462.5 ± 2Atli 


2469.1 ±3.7tJ:^ 


2462.6 ±0.7 


T{Df), MeV 


46.6 ±8.1^1^ 


43 fixed 


49.0 ± 1.4 


h{Df) 


— 1 fixed 


— 1 fixed 




Di(2430)7D? 


1.0 fixed 


1.0 fixed 




D*{2A00)yDf 


1.1 ± 1.1 


1.7 fixed 




Feed-downs/£)2° 


0.3 ±0.4 







Table 2: Results of the simultaneous fit for the yields (N), masses (M), widths 
(V) and helicity parameters (h) of the D\ and D'^ mesons, for the ratios of the 
wide states Di(2430)° and Do(2400)° to the narrow states and D^^ , and for the 
ratio of the feed-down (see text) to the D"^ — )■ D^n^ . The first uncertainties are 
statistical and the second are systematic. The results (HERA II) are compared to 
earlier ZEUS results at HERA I [1] and to the PDG [16]. 
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HERA II 


PDG 


N{Dt D*^n+) 


759 ± 183 




N{D*-^ D*°TT+) 


634 ± 223 






737 ± 164 




M{Dt), MeV 


2421.9 ±4.71?;^ 


2423.4 ±3.1 


r(D+), MeV 


25 fixed 


25 ±6 


h{Dt) 


3.0 fixed 




MiO;-^), MeV 


2460.6 ± 4.4+1^:^ 


2464.4 ± 1.9 


T{D*+), MeV 


37 fixed 


37 ±6 




-1.0 fixed 





Table 3: Results of the fit for the yields (N), masses (M), widths (V) and helicity 
parameters (h) of the Df and Z^2^ mesons. The first uncertainties are statistical 
and the second are systematic. The results are compared to those of the PDG [16]. 





total 


5l 




^3 




M{D^), MeV 


+0.4 
-1.0 


+0.4 
-0.3 


+0.0 
-0.8 


+0.1 
-0.5 


+0.1 
-0.1 


M{Df), MeV 


+1.3 
-1.1 


+0.9 
-0.9 


+0.9 
-0.5 


+0.2 
-0.2 


+0.0 
-0.1 


r(L)°), MeV 


+1.9 
-5.4 


+1.6 
-2.3 


+0.0 
-1.6 


+1.0 
-4.5 


+0.0 
-0.0 


T{Df), MeV 


+5.9 
-3.8 


+4.0 
-3.5 


+0.1 
-0.2 


+4.3 
-1.7 


+0.0 
-0.0 


h{D\) 


+4.6 
-1.8 


+3.1 
-1.3 


+2.4 
-0.3 


+2.3 
-1.3 


+0.1 
-0.1 



Table 4: Total and 81-84^ (see text) systematic uncertainties for the mass, width 
and helicity parameters of the neutral excited charm mesons. 





total 


Si 


62 


S3 


54 


^5 


M{Dt), MeV 


+3.4 
-1.2 


+3.2 
-0.1 


+0.0 
-0.7 


+0.6 
-0.1 


+0.1 
-0.1 


+0.6 
-0.9 


M(L);+), MeV 


+3.7 
-0.8 


+1.7 
-0.5 


+3.1 
-0.0 


+0.4 
-0.2 


+0.1 
-0.1 


+0.9 
-0.6 



Table 5: Total and 61-6^ (see text) systematic uncertainties for the mass, width 
and helicity parameters of the charged excited charm mesons. 
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total, % 




62,% 


^3,% 


64,% 


^6,% 


^7,% 


X'cxtr 

Dl^D*+-K- /D''+ 


+19.2 
-14.5 


+ 16.4 
-12.2 


+6.7 
-0.0 


+3.4 
-7.5 


+0.3 
-0.0 


+1.5 
-2.0 


+6.5 
-0.0 


X'cxtr 


+13.5 
-18.2 


+ 11.9 
-12.9 


+3.7 
-5.0 


+1.2 
-11.8 


+4.9 
-0.0 


+0.9 
-1.5 


+0.1 
-0.0 


X'cxtr 

Df^D+TT~ /D+ 


+25.2 
-17.3 


+ 18.6 
-7.8 


+11.9 
-0.0 


+5.4 
-15.4 


+1.0 
-0.0 


+0.5 
-0.8 


+10.7 
-0.0 




+20.1 
-19.5 


+9.9 
-13.5 


+0.0 
-4.7 


+9.6 
-3.3 


+0.0 
-0.7 


+2.3 
-2.5 


+14.4 
-12.7 


f{c^Dl) 


+15.8 
-18.6 


+ 11.9 
-12.9 


+3.7 
-5.0 


+1.2 
-11.8 


+4.9 
-0.0 


+0.9 
-1.5 


+8.1 
-3.6 


f{c ^ Df) 


+22.4 
-15.1 


+ 16.1 

-9.1 


+8.9 

-0.0 


+4.0 

-10.7 


+0.6 
-0.0 


+0.6 

-1.0 


+12.2 

-5.3 



Table 6: Total and 6i-6j (see text) systematic uncertainties for extrapolated 
fractions, for ratios of the dominant branching fractions and for fragmentation 
fractions of the D\ and Df mesons. 





total, % 


5l,% 


52,% 


53,% 


(54,% 


55,% 


^6,% 


^7,% 


X'cxtr 

D+^D*On+/DO 


+42.6 
-6.1 


+30.5 
-0.0 


+18.3 
-0.0 


+3.7 
-2.6 


+0.0 
-0.0 


+22.2 
-0.0 


+1.8 
-5.2 


+6.0 
-1.9 


X'cxtr 

D*+^D*OTr+/DO 


+24.6 

-14.8 


+ 14.7 

-1.3 


+6.3 
-2.4 


+1.2 
-7.9 


+0.0 
-0.0 


+ 13.5 
-4.6 


+3.5 
-4.0 


+12.5 
-10.5 


•xextr 

D*+^D^>n+/DO 


+18.0 
-8.0 


+ 13.4 
-0.8 


+5.6 
-4.3 


+0.2 
-5.2 


+0.0 
-0.0 


+3.6 
-0.0 


+1.6 
-1.4 


+9.8 
-3.9 




+23.8 

-19.1 


+ 10.5 

-8.5 


+8.3 
-10.0 


+7.0 
-4.7 


+0.0 
-0.0 


+6.9 
-9.1 


+2.7 
-1.9 


+16.9 

-9.3 


fic Dt) 


+42.7 
-7.3 


+.30.5 
-0.0 


+18.3 
-0.0 


+3.7 
-2.6 


+0.0 
-0.0 


+22.2 
-0.0 


+1.8 
-5.2 


+7.1 
-4.4 


fic^n;^) 


+16.7 
-7.1 


+ 12.0 
-0.0 


+1.8 
-0.0 


+0.5 
-5.4 


+0.0 
-0.0 


+8.2 
-1.2 


+2.5 
-2.7 


+7.7 
-3.6 



Table 7: Total and 6i-6j (see text) systematic uncertainties for extrapolated 
fractions, for ratios of the dominant branching fractions and for fragmentation 
fractions of the Df and Da"*" 'me.sons. 
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ZEUS 



ZEUS 




0.14 0.145 0.15 0.155 0.14 0.145 0.15 0.155 

M(Kjrjr,)-M(Kjr) (GeV) M(Kjtirirjr,)-M(KjtJtJt) (GeV) 



Figure 1: The distribution of the mass difference (dots), (a) AM = M^KtttTs) — 
M{Kt[) and (h) AM = MlKirirTrTTs) — M^Kimir). The solid curves are fits to 
the sum of a modified Gaussian function and a background function (dashed lines). 
Candidates from the shaded area, 0.144 — 0.147 GeV, are used for the analysis of 
excited charm mesons. 
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ZEUS 



ZEUS 




1500 



1000 



500 



N(D*) = 39283 ± 452 - 
S(D*)>3 

J I \ \ \ I I \ L 



> 

" 40000 



g, 35000 

S 
_o 

"S 30000 

s 



20000 



15000 



10000 



5000 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


- b) 


• ZEUS (373pb-^) 




— Gauss'""'^+bg. -_ 




— background 












N(D'') = 145740 ±2944 






S(D")>0 


1 1 1 1 1 1 




1 1 1 1 1 1 1 1 1 



1.8 1.85 1.9 1.95 

M(KjiJi)(GeV) 



1.8 1.85 1.9 1.95 

M(Kji) (GeV) 



Figure 2: The mass distributions (dots), (a) M^K'tt^tt^) for events with sig- 
nificance S > 3 and (h) M{K~7r^) for events with significance S > 0. The solid 
curves are fits to the sum of a modified Gaussian and a background function ( dashed 
lines) and for (h) including also a contribution from a second broad modified Gaus- 
sian representing a reflection (see text). Candidates from the shaded areas, (a) 
1.85 — 1.89 GeV and (b) 1.845 — 1.885 GeV, are used for the analysis of excited 
charm mesons. 
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ZEUS 



700 
600 ^ 
500 ^ 
400 ^ 
300 ^ 
200 ^ 
100 ^ 



220 



— 1 — I — I — I — I — I — I — I — I — r- 



a) 



• ZEUS (373pb"') _z 




N(D^")= 1798 ±293 / ', , 



2.1 



2.2 



2.3 



2.4 2.5 2.6 2.7 2.8 

lext 



2.9 



M(D \) = AM'^^+MCD (GeV) 

PDG 



-I— p-i— I— I— I— p-i— I— I— I— p-i— I— I— I— p-i— I— I— I— 
— simultaneous fit 




60 
40 
20 




•■■ signals 

■■■ feed-downs 

- - background 

backgr.&wide-state 



2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 

MiDX) = AM'"'VM(D^pj^^ (GeV) 



Figure 3: The mass distributions (dots), a) M{D*'^TTa) and b) M{D'^TTa)- The 
solid curves are the result of a simultaneous fit to a) and and to b) and 
feed- downs plus background function (dashed curves). The contributions of the wide 
states Di(2430)° and Dq(2400)° are given between the dashed and dotted curves. 
The lowest curves are the contributions of the D\, and feed-downs to the fit. 
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ZEUS 



> 

^ 1200 



200 



N(D*+) = 737±164 



N(D*+) =634 ±223 

2 FD 

N(D^) =759 ±183 

1 FD 



S(D^>0 



• ZEUS (373pb-^) 
fit 

D*2^ signal 

D| feed-down 

feed-down 

background 




Qi...L-.i^.'.^-j-.-i-.l^ki'i'mitl.L2c'jtjJj,**lt,i'iu^-.L...^.^t.Tis.tiM.^^ 
2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 

M(D\)=AM''"+M(D") (GeV) 

PDG 



Figure 4: The mass distribution (dots), M{D^Tia)- The solid curve is the result 
of a simultaneous fit to the feed-down (FD) Df and contributions and to the 
02'^ signal plus background function (dashed curves). The lowest curves are the 
contributions of the Df and to the fit. 
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ZEUS 




2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.2 2.3 2.4 2.5 2.6 2.7 2.8 

M(D*''7iJ=AM""+M(D*"') (GeV) 



Figure 5: The mass distributions (dots), M{D*~^7ia) in four helicity intervals: 
(a) |cosa| < 0.25; (b) 0.25 < |cosa| < 0.50; (c) 0.50 < |cosa| < 0.75; (d) 
|cosa| > 0.75. The solid curves are the result of the simultaneous fit to D\ and 
plus background function (dashed curves). 
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ZEUS 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



rs/(rs+r„) 

Figure 6: The allowed region of cos (f), where is the relative phase of S- and D- 
wave amplitudes, versus the fraction of S -wave in the — t- D*7r decay for ZEUS, 
BABAR and CLEO measurements. 
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